Study of the response of '°B-doped MCP to wide-energy range neutrons from eV to MeV* 
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Neutron-sensitive microchannel plates (nMCPs) have applications in neutron detection, including energy 
spectrum measurements, neutron-induced cross-sections, and neutron imaging. ‘°B-doped MCPs (B-MCPs) 
have attracted significant attention owing to their potential for exhibiting a high neutron detection efficiency 
over a large neutron energy range. Good spatial and temporal resolutions are useful for neutron energy-resolved 
imaging. However, their practical applications still face many technical challenges. In this study, a B-MCP 
with 10 mol% '°B was tested for its response to wide-energy neutrons from eV to MeV at the Back-n white 
neutron source at the China Spallation Neutron Source. The neutron detection efficiency was calibrated at 1 eV, 
which is approximately 300 times that of an ordinary MCP and indicates the success of '°B doping. The factors 
that caused the reduction in the detection efficiency were simulated and discussed. The neutron energy spectrum 
obtained using B-MCP was compared with that obtained by other measurement methods, and showed very good 
consistency for neutron energies below tens of keV. The response is more complicated at higher neutron energy, 
at which point the elastic and non-elastic reactions of all nuclides of B-MCP gradually become dominant. This is 
beneficial for the detection of neutrons, as it compensates for the detection efficiency of B-MCP for high-energy 
neutrons. 
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I. INTRODUCTION imaging with epithermal fast neutrons, is a powerful method 
for studying the microstructures of materials and has become 
Neutron energy-resolved imaging (NERD), such as Bragg- an attractive method in the field of neutron imaging in re- 


edge imaging with thermally cold neutrons and resonance Cent years. For decades, detector technologies for neutron 
imaging, such as scintillator arrays, position-sensitive gas de- 


tectors, as well as CCD and CMOS cameras, have been lim- 

ited by the neutron detection efficiency and could not fully 
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the capability of the NERI method and are probably the most 
cutting-edge method in the field of neutron photography [8— 
11]. 

In 1990, Fraser and Pearson first proposed the nMCP and 
doped °Li,O into MCP to obtain a thermal neutron detection 
efficiency of only 0.21% [12]. Since then, !°B and "Gd 
with larger thermal neutron cross sections have been used 
as the main doping materials [13, 14]. For example, Ts- 
inghua University developed an event-counting thermal neu- 
tron imaging detector using an nMCP doped with 3 mol% 
natGd.O3; the thermal neutron detection efficiency reached 
33% [15]. Harbin Engineering University used mixed doping 
of 10 mol% '°B2O3 and 0.4 mol% "*'Gd2Oz to achieve a 
thermal neutron detection efficiency of 28%, which was better 
than that doped with !°B only [16]. The products a and ‘Li of 
the (n, 1?B) reaction have stopping ranges of approximately 
2-3.5 um in the MCP [17], whereas the product y of the (n, 
nat Gd) reaction has a range of approximately 100 um. There- 
fore, !°B-doped nMCPs have a relatively higher spatial reso- 
lution than those doped with "“*Gd. On the other hand, !°B 
has higher and smoother reaction cross-sections for a neutron 
energy above the epithermal region, which is more suitable 
than "Gd, especially for epithermal-fast neutron resonance 
imaging. Increasing the doping concentration to BOdetection 
efficiency of up to 60% for thermal neutrons. !°,3However, 
until now, the development of MCPs with high concentrations 
of 10B has been difficult because excessive !°B corrodes the 
MCP glass substrate during the process [19]. !°B-doped nM- 
CPs (B-MCP), developed by U. C. Berkeley, NOVA Corpo- 
ration, Tsinghua University, etc., are mainly used in thermal- 
cold neutron imaging [20-22]. However, neutron resonance 
imaging in a wide energy region from eV to MeV has higher 
detection efficiency requirements and must solve the complex 
response problem of B-MCP to high-energy neutrons. 

This paper presents a study on the detection efficiency and 
responses of a B-MCP to wide-energy neutrons from eV to 
MeV at the Back-n white neutron source of the China Spal- 
lation Neutron Source (CSNS). Back-n provides a white neu- 
tron beam with a wide energy spectrum, high neutron flux, 
and good temporal resolution for serving as a multipurpose 
platform to support nuclear data measurements, detector cali- 
brations, and neutron photography [23—26]. The main factors 
affecting the detection efficiency loss of the B-MCP, signal 
sources in the B-MCP output, and the contribution of differ- 
ent reactions to the neutron detection efficiency in different 
energy regions were analyzed. The results provide a basis 
for the ongoing development of B-MCP for neutron imaging, 
especially for energy-resolved imaging over a wide energy 
range. 


Il. EXPERIMENTAL SETUP 


The B-MCP used in the experiment was developed in co- 
operation with North Night Vision, and its main parameters 
are listed in Table 1. 

The B-MCP contains approximately 10 mol% !°B from 
the glass substrate doping, and composition measurements 


Table 1. The main parameters of the B-MCP and an ordinary MCP. 


Parameters B-MCP Ordinary MCP 
Diameter(mm) 33 33 
Thickness(mm) 0.515 0.62 
Body resistance(MQ) 81 146 
Pore inner diameter(um) 8 10 
Pore area ratio (%) 60 60 
Gain 1915@800V__1.3x 10° @1200V 
10B (mol%) 10 = 


show that the preparation process does not induce a sig- 
nificant boron loss. The ordinary MCP also contains trace 
amounts of approximately 0.15% natural boron (19.78% !°B 
and 80.22% 11B) originating from the preparation process 
(not shown in the table). The gain of the B-MCP is approxi- 
mately two orders of magnitude lower than that of an ordinary 
MCP, which is beneficial for neutron detection. A large gain 
in a B-MCP would lead to an excessive amplitude range of 
neutron signals, which is not conducive to the efficient acqui- 
sition of neutron events with specifically designed electronics 
that have a limited dynamic range. Therefore, the B-MCP 
should reduce the gain appropriately and work in unison with 
an ordinary MCP to obtain the required total signal gain. To 
improve the secondary electron efficiency and reduce the ion 
feedback, MCP pores have a bevel angle of approximately 
5°~12° and are stacked in the V-shape of the two MCPs. The 
nuclear reactions of !°B with incident neutrons are expressed 


by Eq. (1). 


10B +n >’ Li(1.0MeV) +* He(1.8MeV) 6% 


10B +n =>” Li(0.83MeV) +* He(1.47MeV) +y 94% 
(1) 


The simulations ( Sec. IV) indicate that a large amount of 
a or Li, which can penetrate the substrate, enter one or more 
pores of the B-MCP and excite the primary electrons on the 
inner wall. When a specified electric field is applied to the 
MCPs, the electrons multiply many times in the pores to form 
an electron cloud at the exit surface of the MCP. The electron 
cloud is generally sufficiently strong such that it can be col- 
lected by the collector electrode. The obtained primary signal 
reflects the response of the B-MCP to neutrons. In this exper- 
iment, a single-anode PCB board was used as the B-MCP sig- 
nal collector electrode; the package structure of the B-MCP 
detector is shown in Fig. 1. 

The main structure includes a ceramic cover, cushion, four 
voltage-partition electrodes, a V-shaped B-MCP/MCP stack, 
a ceramic base, and an anode PCB board. The ceramic cover 
and base serve as external structures that provide integrated 
packaging. Four copper ring electrodes are used to apply volt- 
age to both sides of the B-MCP and MCP, and the cushioned 
copper ring protects the B-MCP/MCP stack from being over- 
squeezed. An ordinary MCP/MCP stack with the same pack- 
age structure is used for comparison. The resistor series are 
used to divide the voltage, and the voltage ratio of the B-MCP 
and the MCP is 0.9:1 for the B-MCP detector and 1:1 for the 
MCP detector. The anode is grounded using a 2-MQ resistor. 


Fig. 1. Schematic diagram of the B-MCP mechanical package struc- 
ture. 


The experimental layout of the B-MCP at the Back-n 
white-neutron source is shown in Fig. 2. The B-MCP (or 
MCP) detector is mounted on a bracket that allowed the neu- 
tron beam to pass perpendicularly through the center of the 
MCP. The surface of the B-MCP is located approximately 58 
m from the spallation target. The beam profile at the detector 
is set to ®20 mm, which is smaller than the diameter of the 
B-MCP, thus allowing the neutron beam to pass completely 
through the sensitive area of the detector. The neutron flux 
is approximately 1.87x 10° neutrons/cm?/s, and the energy 
range covers 0.5 eV-200 MeV [27, 28]. The B-MCP detec- 
tor and preamplifier (Mesytec MPR-1) were placed inside the 
vacuum chamber of the LPDA spectrometer [29]. The exper- 
iment used the common electronics and data acquisition sys- 
tems of the Back-n facility [30, 31]. The output signals of the 
preamplifier were digitized by the readout electronics system 
with a 12 bit folding-ADC and a 1 GHz sampling rate [32]. 
The system was triggered by the TO signal from the proton 
beam monitor, and the threshold was set to 14.6 mV to fil- 
ter the electronic noise. The sampling duration was set from 
2 us before the TO trigger to 8 u s after the TO trigger, and 
the waveform of each neutron event was recorded. Because it 
is difficult for MCPs to distinguish between dark counts and 
neutron events, the B-MCP uses a voltage of -1750 V to keep 
the dark counts below 0.1 cps/cm?. Each measurement took 
10800 s to obtain accurate statistics. 


I. EXPERIMENTAL RESULTS AND ANALYSIS 


Fig. 3 shows a typical signal waveform recorded immedi- 
ately after the TO trigger. The data analysis was performed 
using the ROOT program developed by CERN [33]. Owing 
to the long decay time of the preamplifier output, a signal 
pileup may occur within a sampling window of 10 u second. 
To filter out noise and discriminate between effective events, 
we applied smoothing and differential transformation meth- 
ods to the primary waveforms. The waveform after the differ- 
ential transformation is presented, where the baseline of the 
differential waveform recovers to zero immediately after each 


event pulse. This signal caused a sudden falling edge of the 
primary waveform, corresponding to a negative extremum af- 
ter differentiation. By setting an appropriate threshold, these 
extrema can be selected to obtain signals within the sampling 
window. There were significantly more signals close to the 
TO trigger than far from TO because fast neutrons arrived in 
shorter flight times. 

Fig. 4 shows the time-of-flight (TOF) pulse height (PH) 
correlation spectrum of the total B-MCP signals. Two event 
bands around 195 us and 605 us correspond to the double- 
bunch time structure of proton pulses from the CSNS accel- 
erator and represent the y-flash generated in the spallation 
target. The yflash detected by the B-MCP is also used to 
calibrate the TO signal because the flight time of the flash 
from the target to the B-MCP is 195 us. Small signal clus- 
ters before the yflash were considered random noise, whereas 
the signals between the y-flash bands corresponded to high- 
energy neutrons from the first neutron bunch. Based on the 
noise, we set the signal threshold to a five ADC-code to ob- 
tain pure neutron events for the trend and structural analy- 
sis of the neutron energy-count spectrum. The data acquisi- 
tion threshold of 14.6 mV resulted in the loss of some small 
neutron signals, causing a rare-event gap with a PH below 
20 and TOF above approximately 24000 (corresponding to a 
neutron energy of approximately 30 keV). These events were 
sparse in time distribution and were less likely to be passively 
recorded with an adjacent large-signal trigger than neutrons 
in the higher-energy region. 

The kinetic energy of neutrons can be calculated using Eq. 


(2): 


1 (2) 


where v = Lftignt/TOF. The neutron flight length 
L flight Was Corrected to 58.04 m using the resonance peaks in 
epithermal neutron region. The neutron energy-count spectra 
of the B-MCP and MCP detectors are shown in Fig. 5. 

A large difference between the two curves is observed, par- 
ticularly in the range of approximately 0.5 eV-30 keV, where 
the counts of the B-MCP are two orders of magnitude larger 
than those of the MCP, indicating the successful doping of 
10B into the B-MCP. The trends of the two curves were al- 
most constant in this region and gradually approached above 
30 keV. Notably, the counting rate of the MCP increased sig- 
nificantly and was close to that of the B-MCP for neutron 
energies larger than MeV (reaching approximately 80% of 
the B-MCP). Simulations confirmed that above 30 keV, the 
contribution to the counts from the nuclear reactions by the 
other elements in both the B-MCP and ordinary MCP grad- 
ually increases. This result will be presented in Sec. IV. 
The intrinsic microstructure of the B-MCP curve below 30 
keV reproduced the original energy spectrum of the Back-n 
neutron beam, which was confirmed by the retrieved energy- 
count spectrum of the B-MCP using the !°B (n, œ) cross- 
section from the database [34] as shown in Fig. 6a. The two 


Neutron beam 


Fig. 2. Experimental layout of the B-MCP at the Back-n white neutron source. 
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Fig. 3. A typical waveform recorded close the TO trigger and its differential transformation. 


normalized curves were consistent in both profile and struc- 
tures below approximately 30 keV. The resonance peaks of 
some nuclides such as Ta, Mo, Mn, and Fe in the spallation 
target are marked. This confirms that events below 30 keV 
primarily originate from the doped !°B. Similar structures 
are also observed in the MCP spectrum, which are attributed 
to the small amount of natural boron contained in an ordi- 
nary MCP substrate. A similar energy count spectrum of the 
MCP using the !°B (n, a) section is shown in Fig. 6b. The 
two normalized curves were entirely consistent in the profile 
and roughly consistent in the structures below approximately 
30 keV. However, the significant increase in counts above 30 


keV was inconsistent with '°B (n, a) reactions. Additional 
complex responses were produced by the interaction of MCP 
components with higher-energy neutrons. The two curves in 
Fig. 5 have similar trends and structures above MeV, reflect- 
ing neither the Back-n spectrum nor the (n, œ) cross-sections 
of 10B. These curves reflect the cross-sectional structures of 
the responding nuclides in the MCPs. 


The detection efficiency of the B-MCP was calculated us- 
ing Eq. (3) from eV to MeV, and the results are shown in Fig. 
J. 
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Fig. 4. The TOF-PH spectrum of the total B-MCP signals. 
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Fig. 5. The neutron energy-counts spectra measured by the B-MCP 
and the ordinary MCP. 
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where ®Sgt is the total neutron flux in the B-MCP, nı is 
the ratio of the measured counts in the energy bin to the total 
measured counts, and 772 is the ratio of neutron flux in the en- 
ergy bin to the total neutron flux. We compare the detection 
efficiency with (n, a) cross section of !°B [34] as a function 
of the neutron energy in same span of orders of magnitude. 
The two curves below 30 keV have the same decreasing trend, 
indicating that this part of the detection efficiency is from the 
(n, a) reaction of 9B. The detection efficiency above 30 keV 
does not decrease with the cross-section of the !°B, which is 
due to additional complex responses from the MCP compo- 
nents, as analyzed in the Fig. 6. 

Taking the detection efficiency of 1 eV neutron as a cri- 
teria, the counts with a PH threshold larger than five ADC- 
code corresponding to a neutron detection efficiency of about 
0.38%. Although this is about 293 times larger than the or- 
dinary MCP, it is far below expectation of 8% calculated us- 
ing the !°B doping concentration. Primary waveform analysis 
suggests that a large number of neutron events with PH less 
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Fig. 6. The neutron energy spectrum measured by (a) the B-MCP 
and (b) the ordinary MCP, and their normalized comparison with the 
back-n spectrum. 
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Fig. 7. Neutron detection efficiency of the B-MCP from eV to MeV 
and its comparison with (n, œ) cross section of UR, 


than five are not well separated from the noise. The efficiency 
loss may be due to several reasons. A major contribution is 
that some !°B reaction products cannot enter the MCP pores, 


as the simulations suggest. There are also smaller neutron 
signals that are excluded by the data acquisition threshold of 
14.6 mV. It is the inherent issue due to the gain of the B-MCP. 
The signal amplitude depends on the position where the pri- 
mary electrons are generated in the B-MCP. The large gain of 
the B-MCP leads to thousand-fold difference in the amplitude 
of the signals, resulting in the inevitable presence of small 
signals that are excluded by the threshold due to mixing with 
noise. Reducing the B-MCP gain effectively compresses the 
amplitude range of neutron signals. Alternatively, increasing 
the ordinary MCP gain to saturate the signals can eliminate 
the amplitude differences. However, it is crucial to consider 
the distinction from noise, dark counts, etc. 


IV. SIMULATION AND COMPARISON WITH THE 
EXPERIMENT 


To further explore the detection efficiency and neutron 
response mechanism of B-MCP, the Monte Carlo code 
GEANT4 was employed to simulate the reactions of the B- 
MCP to the Back-n neutron beam. The FTFP_BERT_HP was 
used as the physics list of the process [35]. The components 
of the B-MCP are listed in Table 2. 


Table 2. The mass percentage of the B-MCP components. 


Components Mass(%) 
SiO2 43 
10B203 8 

PbO 24 
Bi203 13 
Al203 1.5 

K20 8 

Na2O 1.5 

BaO 1 


SiO2 is the primary substrate of the B-MCP, comprising 
approximately 43% of its mass. BigO3 and PbO are used 
to adjust the resistance of the MCP to form an electron ex- 
citation layer, while Al2O3 is to eliminate possible bubbles 
in the MCP and improves electron excitation. Other metallic 
oxide components can increase the hardness and strength of 
the MCP glass. The geometry of the B-MCP is illustrated in 
Fig. 8. 

The arrangement of the B-MCP’s pores is hexagonal struc- 
ture similar to honeycomb [36], with an average distance be- 
tween the pores of about 2-3 um, based on the pore diameter 
and area ratio. The actual B-MCP contains millions of pores, 
making it unrealistic to model the entire structure exactly. A 
local model with a diameter of approximately 100 um and 
hundreds of pores is sufficient to meet the simulation require- 
ments. The Back-n energy spectrum [37] is used, and the 
beam diameter is set to 45 um to ensure that all neutrons pass 
perpendicularly through the B-MCP surface. The pores have 
a bevel angle of 8° to conform to the B-MCP object. 

The total responses of the B-MCP are simulated with nat- 
ural abundance of all the elements other than !°B. These re- 
sponses are divided by the (n, tot) cross-sections of the cor- 


Fig. 8. Schematic diagram of the B-MCP geometry. 
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Fig. 9. The neutron energy spectrum simulated by using total re- 
sponses of the B-MCP, and its normalized comparison with the back- 
n spectrum. 
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Fig. 10. The neutron energy-counts spectra simulated with total ‘°B 
(n, a) reactions and with the signal cases in which charged products 
“Li or a entering the MCP pores. 
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Fig. 11. The neutron energy-counts spectra simulated with: (a) total 
elastic scattering responses and their signal cases; (b) total nonelastic 
reactions and their signal cases. 
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Fig. 12. The simulated total signal responses of the B-MCP and their 
normalized comparison with the experiment. 


the neutron energy spectrum detected by B-MCP, which is 
compared with the original Back-n spectrum, as shown in 
Fig. 9 [34]. The two spectra exhibit consistent profile and 


fine-structures, indicating that the simulation results are reli- 
able. However, not all the responses can generate signals in 
the experiment. It is considered that only the charged prod- 
ucts that enter the pores have good opportunity to excite pri- 
mary electrons and produce signals. An experimental event 
requires that at least one or more of the charged products from 
the reaction pass through the inner wall into the pores. 


The main reactions between neutrons and the B-MCP are 
categorized into three types: the (n, œ) reaction with !°B, 
elastic scattering with all the nuclides, and nonelastic reac- 
tions with all the nuclides other than !°B. The neutron energy- 
counts spectra for the total !°B (n, œ) reactions and for the 
signal cases in which charged products ‘Li or a entering the 
pores are shown in Fig. 10. The proportion of signal cases is 
approximately 61% of the total !°B (n, a) reactions, which is 
considered one of the factors contributing to the loss of neu- 
tron detection efficiency in the experiment. It is observed that 
the proportion of signal cases in the energy range of about 0.5 
eV - 300 keV remains relatively constant, which can be at- 
tributed to the fact that product energy is almost independent 
of neutron energy in the low-energy region. However, the 
proportion significantly increases at higher energy, because 
the product particles carry also the kinetic energy of the neu- 
trons and increases the chance of entering the pores. 


The neutron energy-counts spectra for total elastic scatter- 
ing responses and for their signal cases in which recoil nuclei 
entering the pores are shown in Fig. lla. Elastic scattering 
occurs in all neutron energy regions from eV to hundreds of 
MeV, but only cases with neutron energy above tens of keV 
have a chance of recoil nuclei entering the pores to generate 
signals. Recoil nuclei with low-energy neutrons do not gain 
sufficient kinetic energy. The lighter the nucleus, the easier 
it is to obtain higher recoil kinetic energy. The recoil nu- 
clei entering the pores are mainly 160, followed by 2°Si and 
10B, Moreover, 16O and ?°Si are the main components of B- 
MCP. The total nonelastic reaction responses of the B-MCP 
and their signal cases in which secondary charged particles 
enter the pores are shown in Fig. 11b. Nonelastic reactions 
generally have a high reaction threshold and mainly occur 
with the neutron energy above several hundred keV. 


The summation of all the signal cases described above is 
normalized to the experimental result, as shown in Fig. 12. 
Signal responses below 30 keV are mainly from !°B (n, a) 
reactions. Additional responses by the B-MCP components 
above 30 keV observed in the experiment are mainly at- 
tributed to the elastic scattering of 160, accounting for about 
75% of all the elastic scattering signals. Nonelastic reactions 
occur when neutron energies are greater than 0.5 MeV and 
contributes to the majority of the counts above 10 MeV. They 
are mainly !°O (n, a), ?SSi (n, p), 22K (n, p), etc., as shown in 
Fig. 11b. The nonelastic response signals of the above three 
nuclides as the main factor account for about 82% of the to- 
tal nonelastic response signals. The simulation results are in 
general consistent to the experimental one for neutron energy 
above tens of keV. However, the experimental counts are sig- 
nificantly lower than the simulations in the lower energy re- 
gion. This is consistent with the fact that loss of detection 
efficiency mainly occurs in the low energy region in the ex- 


periment. 


V. CONCLUSIONS AND DISCUSSIONS 
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